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Abstract 

The concept of broken symmetry is used to study bifurcations of equilibria and dy- 
namical instabilities in dynamic model of one-mode laser (nonresonant complex Lorenz 
model) on the basis of modified Hopf theory. It is shown that an invariant set of stationary 
points bifurcates into an invariant torus (doubly-periodic branching solution). Influence 
of the symmetry breaking on stability of branching solutions is investigated as a function 
of detuning. The invariant torus is found to be stable under the detuning exceeds its 
critical value, so that dynamically broken symmetry results in the appearance of low fre- 
quency Goldstone-type mode. If the detuning then goes downward and pumping is kept 
above the threshold, numerical analysis reveals that after a cascade of period-doublings 
the strange Lorenz attractor is formed at small values of detuning. It is found that there 
are three different types of the system behaviour as pumping increases depending on the 
detuning. Quantum counterpart of the complex Lorenz model is discussed. 
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1 INTRODUCTION 

Nonlinear dynamics of laser systems, especially, those that exhibit so-called chaotic behavior has 
been the subject of major interest during the last two decades (see, for example, @|)- In the 
semiclassical approximation, when quantum correlation effects are disregarded, description of 
dynamical instabilities in lasers is shown to be closely related to the theory of nonequilibrium 
phase transitions in dissipative dynamical systems ||. From this standpoint, of particular 
interest are the systems that can be reduced to the well-known models of the dynamical system 
theory. One of the models is the famous Lorenz equations: X = a(— X + Y), Y — rX — Y + 
XZ, Z = —bZ — XY that were originally derived in Q and have been studied intensively from 
the end of the seventies (see, for example, [5-9] and references therein). The above system (in 
what follows it will be referred as the real Lorenz model for X(t), Y(t) and Z(t) are real- valued 
functions) was obtained from a set of hydrodynamic equations in the three-mode approximation 
to describe the convective motion of a layer of fluid that is warmer at the bottom than at the 
top, so that a is the Prandtl number; r (controlling parameter) is the Rayleigh number and 
is proportional to the temperature difference; b depends on geometrical properties of the fluid 
layer. 

In order to clarify the relevance of Lorenz-type models to laser physics let us consider 
the simplest one-mode laser equations taken in the following form ||10|| : 



= — (k + iu) ■ b — ig ■ a 

= -(7 + iu a ) ■ a + ig ■ b ■ S (1) 
= (d Q - S)/T -4g b) 

where b is the dimensionless complex amplitude of the electromagnetic field mode; u 
(k) is the frequency (the relaxation constant) of the mode; a is the dimensionless transition 
dipole matrix element; S is the inversion of the atomic level populations; u a is the frequency 
of the atomic transition; <i is the parameter characterizing the intensity of pumping; g is the 
coupling constant; 7 (T -1 ) is the transverse (longitudinal) relaxation constant. 

Eq.(l) is derived within the framework of semiclassical approach to dymanics of the 
system that constitutes a number of two-level atoms (atomic subsystem) interacting with the 
one mode of electromagnetic field (field subsystem), so that b is the averaged field anihilation 
operator descriptive of the coherent part of radiation. Note that the equations can be obtained 
from Heisenberg operator equations for Dicke-type hamiltonian taken in the rotating wave 
approximation after neglecting of quantum fluctuations and assuming no dependence on spatial 
variables pi ITTf. In addition, the result should be supplemented with the relaxation terms 



as well as the term proportional to the intensity of pumping (see Sec. 5 for more details on 
quantum models behind Eq.(l)). 

After making the substitutions: 

t -> 7 t, b = jX/(2g), a = iS Y/2, S = d + S Z, S = ^n/g 2 

and going over to the representation of interaction the system (1) can be rewritten as a complex 



Lorenz model 12 : 
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X = <r(-(l + zA) • X + Y) 

Y = -(1 - iA) -Y + r ■ X + X ■ Z (2) 
Z = -b ■ Z - U{X ■¥) 

where a = k/j, r = d /S , b = (7T) -1 , A = (uj a — o;)/(« + 7) is the frequency detuning. 
Note that X(t) and Y{t) are complex- valued functions, so that the system (2) consists of five 
real equations. 

In the case of exact resonance when A = 0, though the complex Lorenz model is appeared 
to differ from the real one in some respects, it was shown that basically there is no difference 
between dynamics of the system (2) and one of the real Lorenz model [|T^ . 

In this paper our main purpose is to study detuning induced effects in dynamics of the 
complex Lorenz model on the basis of bifurcation analysis. The latter means that we deal with 
stability and existence of certain branching solutions depending on the intensity of pumping. In 
other words, r is assumed to be a control parameter. By contrast to the real Lorenz model, the 
system (2) has a continuous symmetry group (Lie group of rotations in complex planes), so 
that we approach the problem within the unified concept of symmetry breaking. It should 
be emphasized that the reason behind all qualitatively new effects discussed in the paper is 
precisely the symmetry and A can be served as a quantity to measure influence of dynamically 
broken symmetry on bifurcating solutions and their stability. 

In order to discuss the effects in more specific terms, let us outline some relevant results 
for the real Lorenz model (most of them hold at A = 0). 

As a preliminary we comment on stability and bifurcations of equilibria (steady states). 
The null steady state corresponding to the spontaneous emission regime (no coherent radiation) 
is given byX = Y = Z = and is asymptotically stable at r < Tq = 1. This solution 
loses its stability at r = r and there are two asymptotically stable bifurcating steady states 



X = Y = ±yb(r — 1), Z = r — 1 provided that r > tq. At this stage we have the stationary 
points bifurcated from the origin at r = r$ indicating that the regime of spontaneous emission 
is changed to the initiation of laser generation (convective flow for the problem of [|J). This 
bifurcation produces qualitative changes in phase portrait of the system. The analogy between 
transformation of such kind, that can be regarded as a nonequilibrium phase transition, and 
a second order phase transition in thermodynamic systems led to the synergetic concept of a 
phase transition |0J. 

If a > 6+1, there is another critical value of the control parameter r: r c = — — — — — -, 

(7 — 0—1 

such that the above two stationary points become exponentially unstable at r > r c . Note that 
the linearized operator governing stability of the solutions in question (Liapunov's first theorem) 
has a pair of complex conjugate eigenvalues with negative real parts in the neighborhood of 
r c . These eigenvalues are pure imaginary at r = r c . According to the Hopf bifurcation theory 
0, it follows that a new branching time periodic solution can be expected to appear (Hopf 
bifurcation). Stability of this solution is determined by the Floquet exponents: the solution is 
stable (unstable) if it appears supercritically (subcritically). In the case under consideration 
the bifurcation was found to be subcritical |14j]. It implies that the time periodic solution is 
unstable when the pumping exceeded its critical value given by r c . 
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One of the most striking features of the real Lorenz model is the appearance of so- 
called strange Lorenz attractor instantaneously on passing the control parameter r through its 
critical value r c ("drastic" route to chaos). The term 'strange attractor'(or 'chaotic attractor') is 
commonly used for an attracting set that has a rather complicated structure and each trajectory 
within the attractor is exponentially unstable. There is a number of different quantities to 
measure the complexity (stochasticity) of the attractor structure: capacity (fractal dimension), 
information dimension, Hausdorff and Liapunov (Kaplan- Yorke formula) dimensions, K-entropy 
and so on. 

Detailed description of how the strange Lorenz attractor forms is beyond the scope of 
this paper. In brief, this can be understood as being due to the occurrence of a homoclinic orbit 
in the system: as r passes through the value at which the homoclinic 'explosion' takes place, a 
strange invariant set of trajectories is produced, including an infinite number of periodic orbit 
IE i> f|- Note that, in addition, the real Lorenz model is known to exhibit period-doubling 
[ |15|j , intermittency |16| and hysteresis [|I7]] in various ranges of its parameter space. 

The paper is organised as follows: 

In Sec. 2 it is shown that due to the symmetry the null equilibrium state of the complex 
Lorenz model (spontaneous emission) bifurcate into an invariant set of stationary points at 
r = r = 1 + A 2 . Stability of the equilibrium states is studied as a function of detuning. It 
is found that under a > b + 1 there is a critical value of the control parameter r (pumping 
intensity), r c , such that the states of the invariant set become exponentially unstable at r > r c 
and r c is an increasing function of A 2 . 

In Sec. 3 analytical power series Hopf technique is extended on the system invariant 
under the action of a continuous symmetry group to construct bifurcating solutions and to 
investigate their stability in the vicinity of r = r c at A ^ 0. Due to the symmetry breaking the 
bifurating solution is appeared to be doubly-periodic at A ^ 0. It means that nonzero detuning 
results in the appearance of low-frequency Goldstone-type mode related to the motion along an 
orbit of the group. Moreover, the broken symmetry is found to affect stability of the branching 
solution. 

In Sec. 4 the results of numerical analysis are discussed. Dependence of the relevant 
Floquet exponent on A is calculated. It is obtained that there is a critical detuning, A c , such 
that the bifurcating doubly- periodic solution is stable at |A| > A c . It leads to the formation 
of stable invariant torus. Different routes to chaos depending on the detuning are duscussed. 

Concluding remarks and discussion are given in Sec. 5. Bosonic three-oscillator quantum 
system that gives the complex Lorenz model in the semiclassical approximation (quantum 
counterpart) is formulated. 
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2 SYMMETRY OF THE PROBLEM: STABILITY AND 
BIFURCATIONS OF EQUILIBRIA 

Taking r as a bifurcation (control) parameter, Eq.(2) can be rewritten as an autonoumous 
dynamical system with quadratic nonlinearity: 



x = f(x) = Lx + f 2 (x,x) (3) 



where 



x = (x 1 , x 2 , x 3 , x 4 , x 5 ), X = x 1 + ix 2 , Y = x 3 + ix 4 , Z = x 5 , 
L = Df(0) (L y = ^(0)); £(x,y) = £ — -jj- • x n ■ y m 

UXj n m ux n ux m 

Clearly, the system (2) is invariant under the transformations: 

X -> exp(i0) • X, Y -> exp(i(f)) ■ Y, Z -> Z, 

so that Eq.(3) has a continuous symmetry group G of rotations in the X\ — x 2 and x 3 — x 4 
planes, G ~ SU(1): 

f(T(0)-x)=T(0)-f(x), T(0)GG (4) 
Note that there is an additional discrete symmetry in the complex Lorenz model: 

Y -> x, y -> F, Z -> Z, A -> -A . 

The latter implies no dependence on the sign of detuning. For brevity, in what follows A is 
assumed to be nonnegative. 
Eq.(4) gives 

f(T(0)-x st ) = 0=>\Df(x st )\=0, (5) 

where x st is a noninvariant equilibrium solution to Eq.(3) and its orbit, T((f>) ■ x st , produces 
an invariant set of equilibria. The second equation in (5) is obtained by differentiating the 
first one with respect to 0. Along this line we come to the conclusion that matrix of the 
first approximation, that govern linearized stability of x st , is degenerate and its null vector is 

e = K ■ x st , where K = — T(0)L =o is the generator of G. 

acp 

It can be readily shown that the null stationary point related to the regime of spontaneous 
emission loses its stability after the pumping intensity exceeded its critical value 

r > r = 1 + A 2 (6) 

and the solution bifurcated from the origin can be taken in the form 



x 



f = x f = q = yjb- (r- 1 - A 2 ), x$ = 0, xf = A • q, x 5 = 1 + A 2 - r. (7) 

5 



According to Eq.(5), the stationary point (7) provides an invariant set of the steady states 
corresponding to the appearance of coherent light emission (laser generation). In the first 
approximation stability of the state (7) is determined by eigenvalues of the following matrix: 



£>f(x 



st) 



( 


—a 


A ■ a 


a 












-A -a 


—a 





a 









1 + A 2 





-1 


-A 
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1 + A 2 


A 


-1 
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<1 


-A • q 


-q 





-b 


) 



(8) 



As it can be seen from Eq.(5), \L C \ =0 and the kernel of L c is defined by the null vector 
e = (0, 1, —A, 1, 0). Strictly speaking, we cannot conclude on stability of stationary point in 
the first approximation for one of the eigenvalues is zero. So the underlying arguments need 
to be modified. Obviously, since there is another steady state in any neighborhood of a given 
steady state on the orbit, it cannot enjoy the property of being asymptotically stable. Suppose 
that all other four eigenvalues have negative real parts and sketch a way how stability of the 
invariant set can be studied on the basis of symmetry arguments. First, the set is asymptotically 
stable if 



||x(t) - T(^(t))x st || -»• , oo (9) 

where x(t) is a solution of Eq.(3) with initial conditions taken in the vicinity of a sta- 
tionary point; ip(t) is determined from the condition 

||x(t) - T(il)(t))x st \\ -> min 

at given t, so that the left hand side of Eq.(9) defines the distance from x(i) to the invariant 
set. The latter can be written in the form: 



z(t) =T(-^(t))x(t) ±e (10) 

From Eqs(3, 10) it is not difficult to obtain the equation for z(t) and to derive the expres- 
sion for i/j(t). Then the motion along the orbit can be eliminated from the consideration and 
subsequent analysis of the modified system can be made in the first approximation. Following 
this line we arrive at the conclusion that the invariant set is asymptotically stable in the case 
under investigation. Some rather straightforward algebra on the subject was made in for 
A 0. 

If a > b + 1 (this condition is found to be independent of A), there is a critical value of 
the bifurcation parameter, such that the solution in question is unstable at r > r c and L c has 
a pair of complex conjugate imaginary eigenvalues Ai = A 2 = ia at r = r c . Equations for r c 
and a can be derived by making use of Routh-Hurwitz criteria: 



P1P2P3 = p\+Popl (11) 
a 2 = ^, (12) 

P3 
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where pi are coefficients of the characteristic polynomial: 



Po = 2ab(l + a)(r - 1 - A 2 ) (13) 

pi = 6((1 + A 2 )(j 2 + 3a(r - 2) + r) (14) 

p 2 = (a + l) 2 + l + r + 2a + A 2 ((a-l) 2 -b) (15) 

p 3 = 2a + b + 2. (16) 

Dependence of the ratio r c and r on detuning is presented in Fig. 1 at b = 1 for cr = 3 
(solid line) and o = 5 (dashed line). It is seen that r c is an increasing function of A. 

In addition, it can be obtained that 9ft(— Ai| r=r J > and then, if not symmetry induced 
degeneracy of L c , we could meet all conditions for the occurence of the Hopf bifurcation. 



3 HOPF-TYPE BIFURCATION 

There are different techniques to investigate Hopf bifurcations and stability of bifurcating time 
periodic solutions @. In this section we adopt the perturbative method, closely related to 
the analytical technique by Hopf, to explicitly construct the bifurcating solution in the form of 
power series over small parameter e characterizing amplitude of the solution in the neighborhood 
of the bifurcation point r = r c [[/], [19] . Coefficients of the power series can be derived by making 



use of Fredholm alternative and linearized stability of the solution, determined by the Floquet 
exponents, can be studied on the basis of the factorization theorem [M. 

In trying to make analysis along the above line we need to modify the method to bypass 
mathematical difficulties arising from the degeneracy of L c . To avoid equations that have no 
solutions it is assumed that the symmetry is 'spontaneously' broken and the bifurcating solution 
is taken in the form: 

x = T(eO) ■ u(t) = T(eO) ■ (x si + e • z(t)) (17) 
Substitution of Eq.(17) in Eq.(3) gives the equation for u(t) : 

u + efTu = f(u) (18) 

Let us introduce the renormalized frequency fl(e) = a/(l + 7~(e)), so that z(s) = z(Qt) 
is a 27r-periodic vector-valued function. The equation for z(s) reads 

a- (z + 0- (e + e-K-z)) = (1 + r) • (L c z + ef 2 (z, z)) (19) 

where L c is defined by Eq.(8) and the dot stands for the derivative with respect to s. 
The 27r-periodic solution of Eq.(19) can be found in the form of power series in e: 

oo oo oo 

z(s) = £ z n {s) e n ,r=Y, r n e n , 6 = ]T cu n e n (20) 

n=0 n=l n=0 
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As it is seen from Eqs.(7, 8), the quantity q = • (r — 1 — A 2 ) can be conveniently 
chosen as a bifurcation parameter, so that 



q = g c + Ag = g c + ]T g n e n (21) 

71=1 

L c = L^+L^-Aq (22) 
e = K* st = e 0) + • Ag (23) 

where are L c , e at q = q c . Recall that has a pair of pure imaginary complex 

conjugate eigenvalues: L^C\ = icxei, L(°)e 2 = -ia e 2 . (For brevity, the superscript (0) will 
be omitted from the eigenvector notations.) Note that eigenvectors {ej, ... e^} of the operator 
adjoint to and {e , ... e 5 } are mutually orthogonal: (e*, e^) = 5ij. 

According to the standard perturbative technique, insertion of Eqs. (20-23) into Eq.(19) 
yields equations to be solved for z n (s) successively. In the case of n > 0, the equations can be 
written in the form: 

j 2n-2 

{« • J" " 4 0) } *n(s) = £ n) (s) = t ] + E lt ] ■ ™viks + c.c. ] (24) 
ds k=i 

Parameters q n , u n and r n can be found by making use of Fredholm alternative. The 
latter states that Eq.(24) has 27r-periodic solutions if and only if 

<eU n) > = (25) 

(el,f1 n) ) = (26) 

Note that Eq.(25) is the solvability condition for time independent part of z n (s) and Eq.(26) 
is to get rid of so-called secular terms. 

In order to reduce remaining arbitrariness in choice of the solutions, z n (s) are subjected 
to the following additional constraints: 



p2tt 



(e*,z n (s)> = (27) 
(el,z n (s))exp(-is)ds = 0, n> (28) 



In the zero-order approximation we have 

{a.^-4°)}zo(s) = -aw eo (29) 

so that 

Zo(s) — A • ei exp (is) + c.c, ojq = (30) 

where A is a complex integration constant that is determined by the initial condition for zo(s) 
and can be eliminated from the consideration by renormalizing the eigenvector ei. 
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After some straightforward calculations the following results can be obtained: 



<?2n+l = T 2n +1 = UJ 2n = (31) 

a-uj, = 2- \A\ 2 - (e*,f 2 ( ei ,e 2 )> (32) 
Zl ( s ) = b + {A 2 ■ bi • exp (2is) + c.c.} (33) 

- g 2 • K(^A C ) = (34) 
where A c is the eigenvalue of L c , such that X c = ia aX q = q c , 

k = 2 • |A| 2 • { 2 • (ej, f 2 (ei, b )> + (el, f 2 (e 2 , bi)> } - a ■ u x ■ (el, Ke t ) (35) 
and the vectors b , bi are solutions of the equations: 

{2ia-4°5}bi = f 2 (ei,e 1 ) (36) 

- L<°> b = f 2 (e x , e 2 ) - (e*, f 2 ( ei , e 2 )> • e . (37) 

At this stage we get the modification of the Hopf theory with allowance for the symmetry 
breaking, so that bifurcating solution is appeared to be doubly periodic. In other words, there 
are two types of frequencies: the basic frequency fl and the Goldstone-type low frequency mode 
with e6 ~ A 2 . (The latter can be inferred from Eq.(32).) It is of interest to note that another 
distinctive feature of the above results is the last term of Eq.(35). This term, being geometrical 
in nature, is accounted for the broken symmetry and will be shown to be of importance to 
stability analysis. 

In the Hopf theory conclusion on stability of the bifurcating solution can be drawn 



from Eqs.(34, 35) based on the factorization theorem |7], 19, 2C] that states about the stability 
depending on the sign of dt(k): if dt(k) < 0, then bifurcation is supercritical and the time 
periodic branching solution is stable at q > q c \ if 3?(/c) > the solution appears subcritically. 

Note that the factorization theorem should be extended to system invariant under the 
action of a Lie symmetry group. For the system under consideration the theorem can be 
recovered by making use ansatz that looks like (17) and real part of the relevant Floquet 
exponent is proportional to 9ft(/c). More details on the subject is given in [f2~H]. 

As a result, Eqs.(34, 35) are key equations for making conclusion on stability of the 
bifurcating invariant torus. In particular, it implies that the torus is stable at r > r c under 
ft(Jfc) < 0. 

4 NUMERICAL RESULTS 

In the previous section we have studied how symmetry of the system affects Hopf-type bifur- 
cation at r = r c . Our findings are: 

1. An invariant set of equilibria bifurcates into an invariant torus. In other words, the 
branching solution is time doubly-periodic, so that Goldstone-type low frequency mode 
is found to appear due to the symmetry breaking. 
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2. It is found that the sign of dt(k) with k defined by Eq.(35) determines stability of the 
torus. The last symmetry induced term in Eq.(35) implies that the broken symmetry 
affects stability of the branching solution. 

Note that the frequency of the Goldstone-type mode as well as the last term of Eq.(35) tend 
to zero as A — ► 0. 

In Fig. 2 are shown dependencies of dt(k) on detuning (A) for a = 5 (solid line) and a = 
10 (dashed line) at b = 1. It is seen that in both cases there is a critical detuning, A c , at which 
dt(k) changes its sign, so that bifurcation being subcritical at A < A c becomes supercritical at 
A > A c . The latter means that an invariant set of equilibria Eq.(7), corresponding to the laser 
generation, bifurcates into the stable torus as r passes through r c under A > A c . 

To get some idea of qualitative changes of attractor structure in relation to detuning, 
there are three trajectories in 3D 9ft (X) — 3ft(F) — Z space in Figs. 3-5 presented at o = 5, b = 1 
and r = r c + 0.2 for various values of A. Fig. 3 is clearly revealed the attractor as an invariant 
torus at A = 0.5 > A c « 0.41, whereas we have the well-known Lorenz attractor under A = 
(Fig. 5). As is shown in Fig. 4, the intermediate case of A = 0.1 corresponds to an entangled 
structure which is hard to interpret. 

One of the ways to clarify the point is to look at the relevant Fourier spectra. To this 
end, the Fourier spectra |X(u;)| and |Z(u;)| are calculated at A = 0.5 (Fig. 6) and A = 0.1 
(Fig. 7). Notice that ^(o;)) 2 is proportional to the power spectrum of the electromagnetic 
field. 

The Fourier spectrum |X(u;)| for A = 0.5, depicted in Fig. 3, indicates the high frequency 
peak at uj ~ a = 5.5 and the two intensive low frequency peaks at uj = and uj rs 0.04. Since 
the frequency uj ~ 0.04 does not contribute to the spectrum |Z(u;)|, this peak can be attributed 
to the Goldstone-type mode. So, the numerical results are in agreement with ones obtained 
from the theoretical analysis of Sec. 3. As far as numerical analysis is concerned, it should be 
emphasized that, working with relatively small number of points (less than 20000), we are not 
to present the results of high precision calculations, but our calculations has been made with 
reasonable accuracy for investigation of the theoretical predictions qualitatively. 

Coming back to Fig. 4 and looking at the Fourier spectra in Fig. 7, let us recall that, 
according to the theory of Sec. 3 and Fig. 2, the relevant Floquet exponent is pure imaginary at 
A = A c . So changing A from above A c , where the invariant 2D torus is stable, downward (r is 
fixed) we encounter another bifurcation point at A = A c , and the torus is expected to bifurcate 
into a 3D torus, embedded in the 5D phase space of the complex Lorenz model. Computer 
simulation confirms this conclusion. Comprehensive analysis of this secondary bifurcation is 
beyond the scope of this paper (some results on the subject were obtained in p2|). Further 



decrease of A would result in other bifurcations. Taking into account that the spectra of Fig. 7 
are typical of period doublings, it can be suggested that the chaotic attractor forms at relatively 
small A after a cascade of doublings. The irregular Fourier spectrum at A = 0.05, displayed 
in Fig. 8, clearly indicates chaotic dynamics of the system. 
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5 DISCUSSION AND CONCLUDING REMARKS 

In this paper we have studied some detuning induced effects in dynamical model of the single- 
mode laser. The key point of stability analysis and bifurcation theory, presented in Sec. 2-3, 
is that dynamical symmetry breaking must be taken into consideration. It is shown that the 
symmetry breaking results in formation of an invariant set of equilibria, which is an orbit of 
the stationary point given by Eq.(7), at r = r and leads to the appearance of low frequency 
Goldstone-type mode related to the motion along the orbit at r = r c . Moreover, it is found 
that there is the symmetry induced term in Eq. (35), so that stability of the branching doubly- 
periodic solution (invariant torus) is affected by the symmetry breaking. 

Coming back to dymanics of the laser system (Eq.(l)), let us discuss what are the effects 
that could be observed experimentally. It should be noticed that, typically, it is difficult to 
meet the condition r > r c in a single-mode laser and it was just a few experiments with gas 
lasers, where the threshold of dynamical instability was exceeded |23|, Q . Based on the results 
of the above theory (Sec. 2-3) and numerical analysis (Sec. 4), in experimental setup of pi" 



with homogeneously broadened one- mode CO2 laser one could expect three different types of 
the system behaviour as pumping increases: 

• At sufficiently small detuning, A < A , (A ~ 0.07 at o = 5 and 6 = 1) the chaotic 
attractor forms abruptly during the passage of r through r c ; 

• The system undergoes a cascade of doublings before its transition to chaos at A < A < 
A c ; 

• If A > A c , the system does not reveal chaotic behaviour even if r w 10 — 20 r c . 

This gives an insight into why just an oscillatory instability was observed for single-mode 
operation of the laser in |^3, |24|. The effect can be attributed to the off-center pumping that 



was used to ensure the single-mode operation, so that the detuning A was greater than its 
critical value A c . 

From the other hand, recently the real Lorenz equations has been employed to study 



three parameter kinetics of a phase transition [p5j . The model was found to represent the main 



features of a second order phase transition in the case of real order parameter. It seems to be 



straightforward to extend the arguments given in |25j to the case of complex order parameter, so 
that the complex Lorenz model could play an important part in investigation of the kinetics of 
a nonequilibrium second order phase transition. Notice that, according to synergetic approach 
[ |i~3|j , a phase transition is realized as a result of mutual coordination between the complex order 
parameter (X), the conjugate field (Y) and the control parameter (Z). So our results can be 
regarded as an extension of the analogy between nonequilibrium phase transitions and phase 
transitions in thermodynamic systems. 

In conclusion, we give some details on a quantum counterpart of the complex Lorenz 
model. As it was mentioned in Sec. 1, in the semiclassical approximation, the well-known 



N-center Dicke hamiltonian 26 : 
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H D = LU-b + b+^--a z + g-(b + a- + ba+) (38) 

N N 

<J ± = Y, a t, " = Y, a r 

r=l r=l 

can provide the Lorenz equations. Recall that equations of motion must be supplemented 
with the relaxation terms as well as the term descriptive of pumping. As far as the problem of 
quantum chaos is concerned, Eq.(38) cannot be considered as an explicit quantum counterpart 
of Eq.(2). Clearly, the reason is that relaxation and pumping do not enter Eq.(38). 

One way to get rid of the above shortcoming is to use two oscillator representation for 
the Pauli operators that enter Eq.(38) p6|, ^7|: cr + — > a 2 ai, o~ — > afa 2 , where af (cij) is the 



bosonic creation (anihilation) operator of the i-th oscillator. The resulting hamiltonian reads: 

H = uj ■ b + b + uji ■ a^ai + u 2 • a\a 2 + g • {b + af a 2 + ba\ a^) (39) 
The next step is to write master equation for density matrix p in the form [RUl : 



- p = i[H, p] + 70 • L b p + 7i • Li p + 7 2 • L 2 p (40) 

where 

Up = {[b + M - [b,pb + ]} + exp(-M ■ {[b,b + p] - [b + ,pb]}, 

L { p = {[af,aip\ - [a h paf]} + exp (-/^u^) • {[a^afp] - [af,pai}}, 

(3 = l/(kBT ), Pi = l/(/c B Tj); k B is the Boltzman constant. It is supposed that the i-th 
oscillator interacts with its thermostat characterized by a temperature of 7] and the thermostats 
are statistically independent (more details on systems of such kind can be found in ER E3). 



Note that the exact solution of the three oscillator model (Eq. (39)) was recently derived by 



making use of the algebraic Bethe ansatz [31 



Assuming that (bafaj) m (b) ■ (a~[aj) and Ti = T 2 (1^ = ji ■ (1 — exp (— fouii))), Eq.(l) 
can be readily derived from Eq. (40). So we have: 

{b, a, S} <-> {(b), (a 2 ai), - a+ai)}, (41) 

uj a = u 2 - ui, k = 70 • (1 - exp (-/3a;)), 7 = 1/T = 2 Ti, (42) 

d /2 = (n) 2 - (71)1, (43) 

where (n)« = (exp (fliUJi) — l) -1 . Note that for d to be positive it is necessary to meet 
the condition: T 2 /T\ > uj 2 juj\. 

Thus, phenomenological parameters are expressed in terms of microscopic quantities by 
Eqs. (41-43) and Eqs. (39-40) yield the explicit quantum counterpart of the complex Lorenz 
model that can be employed to study the problem of quantum chaos |3"2fl . This work is under 
progress. 
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FIGURE CAPTIONS 



Fig. 1 The ratio of the dynamical instability threshold r c to the laser generation threshold 
r = 1 + A 2 as a function of detuning A at b = 1 for a = 3 (solid line) and a = 5 (dashed 
line). It is seen that r c is an increasing function of A. 

Fig. 2 Dependence of $l(k) (see Eq.(35)) on A at b = 1 for a = 5 (solid line) and a — 10 
(dashed line). In both of the cases the coefficient is shown to change its sign at critical 
value of the detuning. 

Fig. 3 Trajectory in 3D $l(X) — 3?(F) — Z subspace at a = 5, b = 1 and r = r c + 0.2 for 
A = 0.5, r c 23.41. The trajectory is shown to wind up the stable invariant torus. 

Fig. 4 Trajectory in 3D 3?(X) - U(Y) - Z subspace at a = 5, b = 1 and r = r c + 0.2 for 
A = 0.1, r c ss 15.28. 

Fig. 5 The strange Lorenz attractor in 3D 3?(X) — $l(Y) — Z subspace at o = 5, b = 1 and 
r = r c + 0.2 for A = 0, r c m 15. The plots in Figs. 3-5 indicate transition of the attractor 
from the invariant torus (Fig. 3) to the Lorenz strange attractor (Fig. 5) as A decreases. 

Fig. 6 Fourier spectra and |Z(c<j)| at o = 5, b = 1 and r = r c + 0.2 for A = 0.5. There 

are three peaks in the spectrum |X(a>)|: uj = 0, uj ~ 0.04 (see the inset in the upper right 
corner), uj a = 5.5. There is no low frequency splitting mode in the spectrum |Z(a>)|, 
so that the second peak is associated with the Goldstone-type mode (see Sec. 3). 

Fig. 7 Fourier spectra and |Z(cu)| at a — 5, b — 1 and r = r c + 0.2 for A = 0.1. Both 

spectra are typical of doublings. 

Fig. 8 The irregular (noisy) Fourier spectrum |X(cj)| at a — 5, b — 1 and r = r c + 0.2 for 
A = 0.05. 
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